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ABSTRACT 

Mergers of spinning black holes can give recoil velocities from gravitational radiation up to several 
thousand km s -1 . A recoiling supermassive black hole in an AGN retains the inner part of its 
accretion disk. Marginally bound material rejoining the disk around the moving black hole releases 
a large amount of energy in shocks in a short time, leading to a flare in thermal soft X-rays with a 
luminosity approaching the Eddington limit. Reprocessing of the X-rays by the infalling material gives 
strong optical and ultraviolet emission lines with a distinctive spectrum. Despite the short lifetime of 
the flare (~ 10 4 yr), as many as 10 2 flares may be in play at the present time in QSOs at redshifts 
~ 1 to 3. These flares provide a means to identify high velocity recoils. 
Subject headings: galaxies: active — quasars: general — black hole physics 



1. INTRODUCTION 

Simulations of binary black hole mergers show large 
recoil velocities ( "kicks" ) of the final merged black hole 
resulting from from anisotropic emiss i on of g ravitational 
radiation [see, e.g., ICampanelli et all (|2007bl ) for a sum- 
mary and references] . The highest kick so far computed is 
2500 km s _1 for spin a* = 0.8, with both black hole spins 
initially perpendicular to the orb ital angular momen- 
tum and anti-aligned to ea ch other (iGonzalez et al.ll2007l ; 
iTichv fe Marronettil l2007h . ICampanelli et all (|2007bf ) 



predict a maximum recoil velocity of 4000 km s for 
equal mass black holes with maximal spin in this configu- 
ration. For a binary supermassive black hole (~ 10 8 Mq ) 
formed during a galactic merger (jBegelman et al.lll980( l. 
the kick may displace the black hole from th e galactic 
nucleus or eject it entirely (jMerritt et al.ll2004l , and refer- 
ences therein) . For a recoil occurring in an active galactic 
nucleus (AGN) with an accretion disk, the inner disk will 
remain bound to the black hole, providing fuel for con - 
tinued AGN activity (jLoebl l2007t IBonning et al J 12007( 1 . 
Such a 'wandering QSO' might be observed as a QSO dis- 
placed from the galactic nucleus or as a QSO with emis- 
sion lines shifted relative to the galacti c velocity. How- 
ever, AGN rarely show displaced nuclei (jLibeskind et all 
|2006|) . and a search for D oppler-shifted broad line regions 
bv lBonning et al.l (120071) failed to detect any strong can- 
didates for recoil. iBogdanovic et ail (|2007l ) argue that 
accretion by the merging black holes will align their spins 
in a way unfavorable for large kicks. 

It is important to explore additional observational sig- 
natures of high velocity recoils that might help to confirm 
actual cases in nature. We note here that the sudden 
change in the black hole velocity leads to an increase in 
the energy of the disk matter that remains bound to the 
black hole, evaluated in the reference frame of the mov- 
ing black hole. This can lead to a brief but powerful 
releas e of energy in soft X-rays and other fo rms of radi- 
ation ((Shields. Bonning. fc Salvianderl 120071 ) . In {J2 we 
discuss the mass and energy involved in the reforming 
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disk and the time scale and luminosity of the resulting 
"recoil flare." In ij3]wc present a simple numerical simu- 
lation of the disk's response to the kick. In S|4]we discuss 
the observational appearance of the recoil flare, includ- 
ing soft X-ray emission from the primary shocks and sec- 
ondary emission from photoionized gas and heated dust. 
In fj5]wc estimate the rate of occurrence of these events 
and the number of recoil flares that may be observable 
at the present time. 

While this man uscript was undergoi ng revision, a pa- 
per appeared by iLippai et al.l (|2008f ) also predicting 
prompt sho cks in an acc r etion disk around a recoiling 
black hole. ILippai et alj (|2008 :) focus on smaller radii 
and shorter time scales, suitable for searching for optical 
counterparts of gravitational wave detections of merging 
supermassive black holes. Here we discuss the outermost 
bound disk with longer timescales with an eye to detect- 
ing recoil flares currently in play in the AGN population. 

2. ENERGETICS 

For a black hole merger taking place in an AGN, 
the accretion disk fueling the AGN will remain bound 
to the recoiling black hole inside the radius i?b = 
(10 1812 cm)M§ wfooo where the orbital velocity equals 
v kick . Here, M 8 = A/ RH /10 8 M W and v mm = ^/lOOO 
km s . For an a disk ((Shakura fc Sunvaevl 119731 : 
iFrank et al.ll2002f) . the disk mass Mb that remains bound 
is 



Mb 



(10 6 - 60 M Q )a:t /5 M 8 1/4 M 7/10 J R^ 4 



(lO 8 -M )a:t /5 M 8 3 / 2 M o 7 / 1 %- 5 o / o 2 



(1) 



(2) 



where Mq is the accretion rate in solar masses per 
year ((Loebl 120071 IBonning et al.ll2007h . Here M refers 
to the normal AGN phase before the tightening bi- 
nary opens a gap in the disk at small radii ~ 10 r , 



where r„ = GMmilc 2 dMi 


osavlievic & Phinnevl 120051: 


MacFadven & Milosavlievid 


20081 lLoebl 12007ft. The re- 



taincd disk mass can fuel QSO activity over a disk con- 
sumption time td on the order of the viscous timescale 
of the outer bound disk where most of the mass re- 
sides. Once the inner gap has refilled, the post-merger 
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accretion rate will resemble the pre-merger rate, giving 

t d « Mb/ Mo « (10 s yrJal^^Mo- 3 / 10 ^ 2 . ^ 

Consider a recoil directed along the rotation axis (z- 
axis) of the accretion disk. Prior to the recoil, the 
material at radius R orbits with the Kcplerian velocity 
^Kcp, specific angular momentum I = i?«Kcp, and energy 
E = — (l/2)«K ep . After the kick, the material still has 
the same angular momentum around the axis of symme- 
try; but it now has a velocity component u z = — Vidck 
along the negative z-axis in the reference frame moving 
with the black hole. (We use u for velocities in the rest 
frame of the post-kick black hole.) This velocity adds in 
quadrature to the orbital velocity = i>Kcp- The total 
energy per unit mass is now E' = E + (l/2)w^ ick . This 
applies at each radius that remains bound to the hole, 
and therefore the total energy of the disk material has 
increased by an amount -Ekick = E' — E = (l/2)MbV 2 ick . 
Because the angular momentum is unchanged, the bound 
material will eventually settle into circular orbits at the 
original radius and binding energy. The excess energy 
that the material has immediately after the kick must 
therefore be dissipated. We argue below that this dissi- 
pation will in large part involve shocks at ~ Wkick, lead- 
ing to an outburst or "recoil flare" with £fl are = ^E-Ekick, 
where £e is of order unity. The energy associated with 
a disk with the mass given by Equation [5] is 

E &axe = (10 5 ™ eT g )aZt /5 M^M^\ w ^%. (3) 

A more rigorous calculation takes account of the slow- 
ing of the black hole as it shares linear momentum with 
the bound disk mass. When the disk has returned to its 
equilibrium state, the disk and black hole share a velocity 
v z = «kick(l + Mx,/ M-qyi) -1 ', and the- energy dissipated is 
£ kick = (l/2)M b «g ick (l + Mb/MsH) -1 - The kinetic en- 
ergy lost by the black hole provides the energy radiated 
and the kinetic energy associated with the motion of the 
disk in the z-dircction. For Mb/MsH < 1 as considered 
here, most of the energy goes to radiation. Additional 
energy will be imparted to disk material outside i?b that 
is drawn toward the moving black hole but remains un- 
bound. We ignore these complications here. 

Following a perfectly axial kick, the particles at any 
radius will follow non-intersecting elliptical orbits in the 
reference frame of the moving black hole, and a number 
of orbits may pass before the excess energy is dissipated. 
For a kick inclined to the disk axis, the available energy 
will be similar (see 33]); but particles initially at radii 
with i^Kep n ot greatly larger than fkick (i-e., R or order 
i?b ) will have a post-kick velocity in the black hole's 
frame that is a strong function of azimuth in the disk. 
These particles will follow orbits with a range of inclina- 
tions and eccentricities, and the disk will be not merely 
perturbed but seriously disrupted. On approximately 
the orbital period, material from different parts of the 
disk will collide and shock at velocities of order t>kick • 
Marginally bound material will lag substantially behind 
the recoiling hole and then fall back into the disk with 
velocity ~ Wkick- More tightly bound material in the in- 
ner disk will have its orbits pert urbed more modestly and 
may shock more gently if at all (jLoebll2007t iLippai et al.l 
120081) . However, for the mass surface density underlying 
Equation [TJ a majority of the bound mass lies outside 
R = 0.5i?b , and we expect that much of the excess en- 



ergy represented by -Eflarc will be dissipated in shocks 
at velocities of order Vkick (see £j3j . The temperature of 
an AGN disk at relevant radii is < 10 4 K, with a sound 
speed c s < 10 km s _1 . Thus, collisions at even a fraction 
of Ukick are highly supersonic. 
We take the time scale for the recoil flare to be 

We = P b £t = (10 3 ' 4 yi)M 8 v^ t , (4) 

where Pb = 27ri?b/wkick is the Kcplerian orbital period at 
i?b • The factors £e and £t express the uncertainty in the 
radiated energy and time scale. Below we estimate £e ~ 
1 and ^ ~ 2 on the basis of simple numerical simulations. 
In §H we argue that the flare may be observable for £e 
as small as 10~ 2 . 
The power associated with this dissipation is about 

Eflare/iflare Or 

L flare = (10 46 - 1 erg s-Vlt/X^^^oVWCt)- 

(5) 

For comparison, the bolometric luminosity of the AGN 
is 

Lboi = eMc 2 / g = (10 45 - 8 erg s" 1 ) M /g (6) 

for efficiency e = 0.1. The factor f g « 10 _1 allows for 
suppression of the ce ntral accretion rate by the inner ga p 
formed by the binary (jMacFadven fc Milosavljevicl2008f ). 
The ratio is 

t It in0.3t-l.r-l -4/5 wl/2 ,V- 3/10 5/2 

Here and below we take £e = 1- The power of the flare 
can substantially exceed that of the AGN. This is true be- 
cause a large part of the disk mass shocks on roughly the 
orbital period, giving a mass shocked per unit time much 
larger than the central accretion rate, which is driven by 
the comparatively long viscous timescale. This disparity 
in time scale offsets the lower energy per unit mass in the 
shocks, compared with the black hole accretion efficiency 
of - 0.1c 2 . 
The post-shock temperature will be 

T sh ock - (10 6 - 9 K)(i> 1000 ) 2 = (0.7 kcVKooo (8) 

(jOsterbrock fc Ferlandl I2006D . The flare lifetime is 
shorter for higher recoil velocities, but the luminosity 
is higher and the X-rays harder. 

The disk mass estimated in Equation [5] is subject to 
a number of uncertainties. We have extrapolated the 
standard a-disk surface density S to radii where the disk 
is too cool for the assumed op acity to apply. However , 
E is a weak function of opacity (jGoodman fc Tanll200l . 
For larger Mbh and lower v kick , the stability limit Mb < 
Mbh becomes important (jLoebll2007l ). If we let i?n m and 
uii m refer to the radius and orbital velocity within which 
the a-disk mass equals Mbh, then from Equation [T] we 
have 

Rvm = (ID 18 ' 11 cra)a 16 { 25 M 3 s /5 M- 14/25 (9) 

and 

wun, = (ID 3 ' 01 km s->lf 2 X /5 Mo /25 . (10) 

This limits the energy and luminosity of the flare, but it 
remains a powerful event. 

A further complicati on is the role o f the disk's self- 
gravity at larger radii (jGood man 2003). For an a-disk, 
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this comes into play where the Toomrc stability parame- 
ter reaches unity, tq fh 10 3 ' 6 r g and a Keplerian velocity 
vq w 10 37 km s _1 . The consequences for the disk struc- 
ture outside tq are unclear. If the surface density of the 
disk is severely reduced at radii ~ 10 5 r g of interest here, 
then the amount of shocked gas and the power in the 
flare will be correspondingly reduced; but clumping of 
the disk gas should not greatly affect our conclusions. If 
the disk regulates itself at Q w 1, then Goodman finds 
a radial dependence of £ oc r _1 or r -3 / 2 rather than 
the r~ 3 / 4 dependence underlying Equation [TJ Such a 
modification of the disk mass near Pb would affect our 
quantitative results but leave the recoil flare a readily 
observable phenomenon. Goodman summarizes a vari- 
ety of observational indications that gas disks in AGN 
do in fact extend to parsec scales. We proceed on the 
assumption that a gas disk exists at these radii, but we 
caution that the mass is uncertain. 

3. NUMERICAL SIMULATION 

In order to test our assumption of high velocity shocks, 
we have a carried out a simple numerical simulation of 
the response of the disk to the black hole recoil. We 
simulated the initial disk as a collection of collisionless 
particles in Keplerian orbit in the x — y plane, and fol- 
lowed the motion of the particles and black hole in 3 
dime nsions with the aid of an N-body code by ([Aarsetbl 
|1985[) . adapted to our purposes. We worked in dimen- 
sionlcss units with G = 1, a particle mass m p = 1, and 
Mbh = 10 5 . The initial disk had zero thickness. (An 
AGN disk at the radii of interest has H/R « 1; and 
for the non- axial kicks of interest, post-kick motions out 
of the original plane far exceed the initial disk thickness 
and internal sound speed.) We ignored the gravitation 
force of the disk on itself and on the black hole. We 
modeled the case of a kick in the y-z plane at 45 degrees 
to the rotation axis. In the units of our model, the ra- 
dial limits of the disk were R\ = 10 and P 2 — 50. We 
took i>kick = 65, giving Pb = 23.6, intermediate between 
Pi and i?2 . Test calculations reproduced a stable disk 
around a stationary hole, and the nonintcrsccting nature 
of orbits for an axial kick. Note that the choice of parti- 
cle mass and radius scale arc arbitrary. The simulation is 
defined by the direction of the kick and the ratios R\ /Pb 
and Rij R\,. Scaling to specific astrophysical parameters 
is straightforward using the expressions in Sj2] 

In order to test for collisions, we associated a fixed 
radius r p of unity with each particle. The initial disk 
comprised concentric rings of particles, separated by 2r p 
in radius and in azimuth. Thus, the disk particles are 
initially just touching their neighbors. This arrangement 
gave us about 2000 particles. We defined a 'collision' to 
occur when two particles passed within r p of each other. 

The expected behavior can be visualized in the rest 
frame of the black hole just after the kick. The par- 
ticles have velocity components — Ukick,^ and — i>kick,z 
added to their original orbital velocity. The total y- 
component in this frame is u y = — VKcp sin(0) — i>kick,2/j 
which varies strongly with azimuth 0. The binding en- 
ergy per unit mass in the frame of the hole is now 
E = —GMbh/R + (1 /2)u 2 , where u is the particle's ve- 
locity in the rest frame of the hole. This determines 
the orbital semi-major axis, a — 2GM/(—E), and period 



P = lis a I (GM)" 1 / 2 . For all radii in our simulation, 
the material in a given ring is bound to the hole for a 
range of <$> and unbound at other azimuths. 

In the units of our model, the orbital period at Pb 
is Pb = 2.28, and this sets the time scale of the recoil 
flare. We ran our simulation from t = to t = 14. Fig- 
ures [TJ - [3] show the positions of the particles at several 
times during the evolution. The more tightly bound ma- 
terial at smaller radii follows the hole relatively closely 
More weakly bound material lags before ultimately re- 
turning to periastron, where it is likely to collide with 
other portions of the disk material. The bound particles 
experience collisions at high velocity with other parti- 
cles originating at substantially different locations in the 
original disk. We take the energy dissipated per collision 
to be (l/2)w 2 ol m p . 

Figure 0] shows the power per unit time produced by 
the colliding particles, computed by summing the energy 
for 100 consecutive collisions and dividing by the elapsed 
time. The power level rises over a time ~ Pb, and per- 
sists for several times Pb . The power at later times is 
exaggerated because of repeat collisions involving par- 
ticles that have completed multiple orbits. The r.m.s. 
relative velocity of the colliding particles, weighted by 
the collision energy, rises from w/i>kick = 0.3 near the 
start of the run to 0.9 at t = Pb and 1.4 for the entire 
run. This characterizes which shock velocities dominate 
the power and the shock temperature and spectrum. 

Two theoretical light curves are shown for comparison 
in Figure |H (1) Here each bound particle was assumed 
to shock exactly one orbital period after the kick, releas- 
ing an energy (l/2)m p v^ ick . The events were ordered in 
terms of P and binned in groups of 50 particles to cal- 
culate the luminosity. (2) The smooth curve in Figure 
H] represents an analytic solution for an axial kick. The 
post-kick orbital period increases monotonically and ap- 
proaches infinity as R approaches Pb ■ We assume that 
the mass originating at each radius shocks at velocity 
Ukick at time t = P, giving a power L = (l/2)v kick dM/dt, 
where dM/dt = dM/dP = 2itY>RdR/dP, and P{R) is 
given by the post-kick binding energy in the frame of the 
moving black hole. For constant surface mass density S, 
this gives 

M = (4/3)(M b /P b )C" 5/3 /(r 2/3 + I) 3 , (11) 

where ( = t/P^. This decreases as M oc i -5 / 3 at late 
times, a qualitatively similar decline to comparison curve 
(1) above. This case is given for reference only, as an ax- 
ial kick will dissipate slowly (see above); and for smaller 
radii (having P << Pb) the orbits are less radically per- 
turbed. Despite the approximations, the various light 
curves agree at the factor-of-two level as to the power 
expected at times of one to several times Pb . 

Integrated over the run, the power dissipated in col- 
lisions is E co \\ = 2.8 x 10 6 , somewhat exaggerated by 
repeat collisions at later times. We may compare this 
with the overall energetics of the disk. Of the 1968 par- 
ticles in the model disk, 831 remain bound to the hole, so 
the bound mass is Mbo = 831. This compares with 381 
particles that are initially inside Pb ; for the tilted kick, 
a substantial part of the mass outside Pb remains bound 
to the black hole. The binding energy of the bound mass 
immediately after the kick is Pbo = —1-28 x 10 6 , and its 
total angular momentum is Z^o = 1-07 x 10 6 . Let us as- 
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sume conservation of angular momentum for the bound 
mass in the aggregate as it settles into a circular disk 
around the recoiling hole. The shallowest binding en- 
ergy for this mass and angular momentum occurs for a 
narrow orbiting ring, giving a radius R IU ig = 16.6 and 
energy E ring = —2.51 x 10 6 . In this case, the energy 
that must be dissipated to reach this configuration is 
AE = E h0 - E ling = 1.23 x 10 6 . Alternatively, if the final 
state is a uniform disk with a sharp outer boundary, then 
Sfmai = (32/25)£ , ring = -3.21 x 10 6 and the dissipated 
energy is AE = 1.93 x 10 6 . These values are comparable 
to the simple estimate AE = (l/2)M b0 Vkick = 1-76 x 10 6 , 
using here the actual bound mass for the tilted kick. This 
supports a value (e~1 in Equation [3] 

For our model, the most tightly bound quarter of the 
particles have a post-kick orbital period P < l.lPb, one 
half of all the bound particles have period P < 2.4Pb, 
and the most weakly bound quarter of the particles have 
P > 5.4Pb- This suggests that, to characterize the main 
surge of the flare, we may take £t ~ 2 in Equation 01 
There is some material with arbitrarily weak binding en- 
ergy after the kick. This material lags far behind the 
moving hole, giving a prolonged tail to the impact rate. 
In our numerical model, the bulk of the shock power cor- 
responds to collisions occurring at radii near or inside i?b 
; the mean collision radius, weighted by collision energy, 
is ~ 12 or 0.5i?b> with little variation during the run. The 
collisions typically involve pairs of bound particles that 
originate at substantially different radii near or outside 
i?b , and substantially different azimuths. These par- 
ticles move on quite eccentric orbits, and the collisions 
producing most of the power typically occur not far from 
periastron. 

We conclude that a non-axial kick will severely disrupt 
the orbits of the mass elements in the disk at radii ~ _R b , 
leading to strong shocks. A physically realistic hydro- 
dynamical simulation will be needed to compute a more 
detailed light curve. However, our simulation supports 
the approximate energetics outlined in &\ 

4. APPEARANCE 

Equation [5] suggests that shocked gas in the post-kick 
disk will release its energy largely in the form of thermal 
soft X-rays, with kT ps 0.7 keV for w kick = 1000 km s -1 , 
or 0.2 keV for ^kick = 500 km s _1 . Subject to the red- 
shift and absorption by the surrounding material, this 
radiation could be observable. 

Consider a QSO with Mbh = 10 8 Mq accreting at 
lO" - 5 Mq yr -1 (giving L/L Ed « lO" 8 ) before the 
inner gap forms. Then vy lm = 750 km s _1 , and for a 
kick velocity of 1000 km s -1 the bound disk mass of 
10 7 ' 7 Mq is less than Mbh- The radius of the retained 
disk in this case is i?b = 10 181 cm, and the surface 
density is S w M h /wRl = 10 42 g cm" 2 . The above 
discussion suggests that, for times t > Pb, the recoil 
flare has the character of material on highly eccentric 
orbits impacting the disk at radii ~ i? b . For the pa- 
rameters just assumed, ifl are « £t-Pb = (10 3 ' 4 yr)£ t or 
~ 5000 yr for £ t = 2. The corresponding mass infall rate 
is M flaro « M b /We = (10 4 - Mq yr" 1 ). This large rate 
of mass undergoing shocks is, of course, the basis of the 
large power of the flare. The density of the infalling ma- 
terial is subject to uncertain geometrical details, but a 



simple estimate is p « Mfl arc / (2ttR^ v) = 10 15 ' 2 g cm 3 , 
or N = pi ran = (10 86 cm -3 ). The cooling time will 
be ~ 10 yr, so that the material will cool in a rel- 
atively narrow zone on the scale of i?b • The local 
energy flux in the X-ray flare is Ffl aro w (l/2),cw 3 s» 
10 8 5 erg cm -2 s _1 . The column density in the impact- 
ing stream is ~ pRb ~ 10 29 g cm -2 . The optical depth 
due to photoelectric absorption by H, He, and heavy el- 
ements will be large, and the X-rays may be unable to 
penetrate the enveloping inflow. However, the luminos- 
ity of the shocked flow, -Lfl are « (10 45 4 erg s -1 ), is close 
to the Eddington limit. This raises the possibility of ra- 
diation driven outflows, although the energy available is 
insufficient to expel the entire disk mass. Given these 
outflows and the complex geometry of the infalling ma- 
terial, there may be clear lines of sight on which the 
X-rays escape. While the shock temperature is high, 
the effective temperature corresponding to Ffl are is only 
T eff = (Eflarc/o-) 174 « 10 32 K. This underscores the 
sensitivity of the observed appearance of the flare to the 
nature of the reprocessing that occurs in the surrounding 
material. 

For an estimate of the spectrum of the cooling shocked 
gas, we computed a coronal equilibrium model using 
version 07.02.00 of the phot oionization code CLO UDY, 
most recently described by iFerland et al.1 (|1998? ). We 
used T = 10 7 K, N = 10 9 6 cm -3 , and solar abundances. 
(This corresponds to a pre-shock velocity 1100 km s -1 
and density of 10 9 cm -3 .) Much of the total cooling is 
thermal brcmsstrahlung. There are many strong emis- 
sion lines, including lines of Fe XVII to Fe XX at 12 
to 15 A, each with several percent of the total energy 
(see Figure [5]). If not totally absorbed by the infalling 
gas, this should be recognizable as a unique spectrum 
for a point source in a galactic nucleus. Comptoniza- 
tion should be insignificant, and the flare will not emit 
hard X-rays. There could, however, be hard X-rays if the 
AGN luminosity is not entirely shut off by the inner gap 
in the disk at the time of the merger. 

The X-ray luminosity will photoionize the infalling 
material approaching the shock. The local energy flux 
in the X-ray flare directed into the infalling gas is 
fflare ~ 10 8 ' 8 erg cm -2 s , giving an ionizing photon 
flux of <pi fa 10 18 ' 2 cm -2 and an ionization parameter 
U = 4>i/Nc = 0.05. We have computed photoionization 
models with CLOUDY for this precursor zone. We used 
the coronal model emission as the ionizing spectrum with 
flux 10 8 8 erg cm -2 s -1 , a gas density N = 10 9 cm -3 , 
and a cutoff column density 10 23 5 cm -2 . The photoion- 
ized depth in this pre-shock flow is ~ 10 13 ' 7 cm, again 
narrow compared to Rb- The high ionization parame- 
ter and hard ionizing spectrum give strong ultraviolet 
emission lines of highly ionized species (see Table Q] and 
Figures [6] and [7]) . The Lya line carries 11% of the in- 
cident energy flux, and O VI and C IV are comparably 
strong. The H/3 and C iv lines respectively have equiva- 
lent widths of ~ 780 and 6100 A in terms of the diffuse 
continuum of the precursor alone. The Balmer and he- 
lium lines are affected by large optical depths. Based on 
the coronal equilibrium model, the shock will contribute 
an additional 12% to the continuum at H/3 and 30% at 
A1550, slightly reducing the equivalent widths. The high 
density suppresses forbidden line emission. Aside from 
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this and transfer effects in some lines, the spectrum is 
not strongly sensitive to density, because the shock ve- 
locity alone determines the post shock temperature and 
also the ionization parameter of the precursor. 

We also computed CLOUDY models for a case corre- 
sponding to «kick = 500 km s -1 . We took M 8 = 10 -05 
and Mq = 10~ 07 to conform to the self-gravity limit. 
This gives ifl are = 10 41 yr for £ t = 2, and a preshock gas 
density of N = 1Q S1 cm -3 estimated as above. Table 1 
shows the emission-line intensities of the precursor ion- 
ized by the emission from a 10 6 ' 3 K coronal model for the 
shocked gas. Many of the same strong emission lines oc- 
cur, such as H/3, Lya, and C iv. However, O vi is much 
weaker, due to a lower level of ionization and a lower 
electron temperature. Thus, the high ionization line in- 
tensities offer a potential diagnostic of the shock velocity 
and hence Ukick • Forbidden lines such as [O in] A5007, 
[O III] A4363, and [Nc v] A3426 have intensities of 0.5 to 
1.0 times /(H/3) because of the lower density. 

We may estimate the relative luminosity of the recoil 
flare at specific wavelengths with the aid of Equation [71 
In the optical, a typical AGN energy distribution gives 
AL A (5100) = O.llLbol (|Kaspi et al J 120001 ). The pho- 
toionized precursor by itself has AL.\(5100) = 0.03Lj^ re 
based on the CLOUDY model. Using these factors in 
equation [7j we have 

Air c /AL AGN = aAC -l / - la= 4/5 M l/ 2 ^-3/10 u 5/2 0) 

(12) 

with C\ = 10 3 for A5100. A similar calculation for 
A1550 gives C\ = 10~ 06 . Therefore, the flare contin- 
uum is comparable with the AGN continuum if f g = 1 
or dominant if / g is small. This in turn implies that the 
equivalent widths of the flare emission lines remain large 
relative to the total observed continuum, including the 
AGN component. In the composite QSO spectrum of 
(jVanden Berk et "all 1200 If ), the equivalent widths of H/3, 
C iv, Lya, and O vi are 46, 24, 93, and 10 A, respectively, 
mostly in the broad line component. With their large 
equivalent widths and narrow profiles, the flare emission 
lines should be prominent in the total spectrum. Indeed, 
they remain conspicuous even if the flare contribution is 
two orders-of-magnitude weaker than estimated here, rel- 
ative to the AGN. The number of detectable flares in play 
could be increased by as much as an order of magnitude 
by the prolonged arrival of weakly bound material catch- 
ing up to the black hole (see Ej3j) . The line widths will be 
~ Wkick, typically 500 to 1000 km s _1 for flares likely to be 
observed (see below). These line widths resemble those 
of Narrow Line Seyfert 1 (NLS1) objects or the lines from 
the traditional narrow line region (NLR) of AGN. How- 
ever, the lines of the recoil flare will be distinctive, show- 
ing (1) potentially large equivalent widths, (2) weakness 
or absence of forbidden lines such as [O in] A5007, (3) 
strong He II A4686 and possibly O VI A1035, and (4) 
likely velocity shifts and asymmetries of order Vkick ■ Ab- 
sorption by dust in the surrounding disk material may be 
important. The disk temperature near i?b is low enough 
before the kick for dust to exist, and the energy flux from 
the shocks may not be enough to evaporate refractory 
grains. 

We may similarly estimate the prominence of the soft 
X-rays from the primary shocked gas. For the AGN 
continuum at 0.3 kcV , Laor et al. (1997) find an av- 



erage vL v q.z = 10 1( V_L„(3000 A); and by the above 
conversions, this gives vL Ut o.3 = 10~ 1S L^ N . From 
the CLOUDY model for T = 10 7 K, we find R . 3 = 
v ^o%/ = 10~ 10 . Combining these expressions 
with Equation we find 

"L%Z/uL$$° = (13) 

io i -Mi?o.3/io- 1 )crV g - 1 a:f 5 M 8 1 / 2 M - 3 / 1 % 1 5 / 2 . 

If a substantial fraction of the X-ray radiation from 
the shocks is able to reach the observer, the object 
will stand out as an exceptionally strong source in soft 
X-rays, relative to its optical and ultraviolet luminos- 
ity. This remains true for fkick = 500 km s _1 , giving 
?shock = 10 6 ' 3 K. Our coronal equilibrium model at this 
temperature shows a continuum dropping at A < 30 A 
and many emission lines forming a pscudocontinuum at 
30 to 60 A. In the vicinity of 40 A (0.3 kcV ), the power 
in the lines corresponds to i?o.3 = 10~ 05 , again giving a 
large ratio of soft X-rays to optical continuum. 

The unified model of AGN posits a thick dusty torus 
that surrounds the emitting core at a radius of one or 
several parsecs (Antonucci & Miller 1985). This torus 
will intercept and reprocess a substantial fraction of the 
X-ray and optical-ultraviolet radiation from the reform- 
ing accretion disk. A large part of the flare luminosity 
should therefore appear as infrared radiation from dust. 

5. RATE OF OCCURRENCE 

The number of recoil flares currently observable de- 
pends on the rate of mergers, the probability of kicks 
of various velocities, and the flare duration at a given 
velocity. Let dN/dt be the merger rate per year of ob- 
server time in some range of Mbh and redshift z\ and 
let f v be the fraction of mergers giving a kick velocity 
greater than v. The number of kicks currently observ- 
able is roughly iV flarc w (dN/dt) [ta aie (l + z)} f v , where 
the factor 1 + z allows for the observed duration of a 
flare with intrinsic duration ifl are . For mergers of two 
black holes with equal spin parameter a* = J/M = 0.9 
and roughly equal masse s satis fying q = TO2/toi>0.2, 
ISchnittman fc Buonaniiol (|2007t ) find / 500 = 0.31 and 
Ziooo = 0.08. These results rely on uncertain assump- 
tions about the merger details but provide an approx- 
imate basis for estimating the number of high velocity 
recoils. 

Tw o alternative e s timates of N n are are as follow s: 

(1) iSesana et all (l200i 120071 ) and iVolonteri I (I2007T ) 
compute binary merger rates in hierarchical merger sim- 
ulations normalized to the observed QSO luminosity 
function. Volontcri (2007, personal communication) has 
kindly provided merger rates for mass ra tios q > 0.2 
for th e various seed black hole scenarios of ISesana et al.l 
(|2007h . In the "VHM" scenario (small seeds), the event 
rate per observer year for 1 < z < 3 is ~ 0.08 yr -1 for 
(mi +TO2) in the range 10 7 to 10 8 Mq and 0.014 yr -1 in 
the range 10 s to 10 9 Mq. The event rate is about 3 times 
smaller in the other scenarios. The rates drop rapidly 
for z < 1. Let /qso be the fraction of these mergers 
that are sufficiently gas-rich to fuel QSO activity. For 
the 10 7 to 10 s Mq range, we have ifl arc « 800£ t yr at 
Wkick = 1000 km s" 1 for M BH = 10 7 5 M . Then in 
the VHM scenario we expect to see roughly 15^/qso 
events currently in play. For i>ki c k > 500 km s , we 
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expect ~ 500£ t /qso events in play, taking account of 
the longer flare duration and the larger fraction of re- 
coils above the lower kick velocity. For 10 s to 10 9 M , 
the limit Mb < Mbh comes into play for a typical value 
L/L Ed = 10~ 05 (e.g., Salviander et al. 2007), and this 
affects ifiaro . If we take tfiare ~ 27ri?ii m /?;ii m (equations 
[5] and [TUJ) , then we expect ~ 8 /qso flares in play above 
1000 km s _1 and ~ 33 /qso above 500 km s _1 . These 
numbers could be several times larger if the disk extends 
beyo nd -Riim with a mass remaining below Mbh- 

(2) iHaehneltl l)1994h derived the black hole merger rate 
from the QSO luminosity function assuming L = L^d, 
one event per black hole, and a QSO lifetime of £qso = 
10 7 ' 6 yr (the "Salpeter time" for black hole growth). At 
Mbh between 10 7 and 10 8 M Q , Haehnelt finds 0.016 
events per year of observer time for z < 3, giving roughly 
3£t flares in play for i>ki c k > 1000 km s" 1 , and 94£ t in 
play above 500 km s~ x . Uncertainties in this estimate 
of iVfl are include the assumed value of L/LEd and the 
QSO lifetime, the flare duration, and the number of ob- 
scured QSOs. However, we do not need the factor /qso 
because the calculation is based on observed QSOs. 

Despite the uncertainties, it appears likely that several 
flares above 10 7 M Q are currently observable. We em- 
phasize that this rests on the assumption that the black 
hole merger occurs while the QSO accretion disk is in 
place. 

For a 10 s M Q hole in a QSO shining at 0.3LEd at 
z = 2 (prior to formation of the inner gap), the re- 
ceived flux from the flare is ~ 10~ 128 erg cm~ 2 s _1 
for £t = 2. The observed spectrum will be softened by 
the redshift, but a considerable fraction of the radiation 
should still be received at photon energies above 0.2 keV. 
At lO" 12 8 erg cm" 2 s" 1 . iHasinger et all ([2005h find - 2 
soft X-ray sources (0.5 to 2 keV ) per square degree, so 
that the recoil flares would need to be identified from 
among roughly 10 5 other sources of similar flux. The 
unique spectrum of the flare should aid recognition. 

An optical search for recoil flares could target the 
strong, narrow emission lines in the rest optical and ul- 
traviolet. For £qso = 10 7 - 6 yr, the fraction of QSOs in 
the recoil flare stage from Equation fj] is 

N Raio /N QSO = 10- 3 - 2 M 8 «r 3 00 £ t /,. (14) 

This gives N &aic /N QSO = 10~ 4 - 3 M 8 £ t for vi oo = 1 and 
10~ Mg£t for wiooo = 0.5. Thus, recoil flares may be 
observable in about one in 10,000 QSOs, and one or 
more examples may be present in existing large spec- 
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troscopic surveys. For examp le, the SDSS DR5 QSO 
Catalog (jSchneider et al.ll2007l ) contains 77,429 objects 
with a median redshift z = 1.48 over ~ 5740 deg 2 of sky. 
About 10 4 of these objects have z < 0.8 so that H/3 is 
accessible, and many have z > 1.7 so that C iv is ac- 
cessible. The higher redshift objects typically have rela- 
tively high luminosity and black hole masses ~ 10 85 Mq 
or higher, so that the self-gravity limit on the accretion 
disk is significant. However, the discussion in $4]suggests 
that, even if the disk mass around i?b is much less than 
assumed in Equation the flare emission lines still could 
be prominent. The observed QSO population serves as 
the parent population for the recoil flares, but the AGN 
luminosity may be dimmed several magnitudes by the 
inner gap at the time of the recoil. This will complicate 
the detection of AGN undergoing recoil flares and their 
inclusion in QSO surveys. 

6. CONCLUSION 

Merged black holes with high recoil velocities can carry 
a large mass of gas in an accretion disk that remains 
bound to the black hole. The kick velocity gives the 
disk material a boost in energy, relative to the moving 
hole, that will in substantial part be dissipated in shocks 
within a few orbital periods as the more weakly bound 
material rejoins the moving disk. This results in a soft 
X-ray flare lasting thousands of years with a luminosity 
rivaling or exceeding the AGN luminosity. The X-rays 
will be heavily absorbed, at least along some lines of 
sight, and the energy rc-cmittcd in line and continuum 
emission from photoionized gas and infrared radiation 
from heated dust. A number of recoil flares may be cur- 
rently observable. Observational detection of such an 
event would provide an important confirmation of the 
occurrence in nature of high velocity recoils. Detailed 
simulations of the observational appearance of these re- 
coil flares are needed to aid in identifying such objects in 
surveys at X-ray and other wavelengths. 
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TABLE 1 

Line Intensities of Shock Precursor 



Ion 


A(A) 


T7.0 


T6.3 


H/3 


4861 


1.000 


1.000 


II(V 


6563 


5.0 


6.5 




1215 


28.7 


88.5 


He I 


5876 


0.10 


0.24 


He II 




17 1 


1 Q Q 

lit . y 


He II 


1640 


2.7 


5.4 


He II 


4686 


0.31 


0.68 


C II 


2326 


1.83 


0.69 


C III 


1909 


5.1 


14.9 


C III 


977 


2.8 


2.9 


C IV 


1549 


25.3 


55.3 


N IV 


1486 


1.58 


5.1 


N V 


1240 


3.52 


1.65 


O III 


1665 


2.0 


4.7 


O IV 


1402 


2.1 


4.8 


O V 


1218 


8.6 


1.6 


O V 


630 


5.1 


0.03 


O VI 


1035 


48.0 


0.44 


Mg II 


2798 


5.2 


5.9 


Si III 


1397 


1.6 


3.6 


log N H 




9.0 


8.0 


log -Fine 




8.77 


6.70 


log F(H/3) 




6.34 


4.21 


EW(H/3) 




780 A 


660 A 


EW(C iv) 




6100 A 


11,500 A 



Note. — Emission-line intensities relative to H/3 of CLOUDY models for shock precursor zone photoionized by coronal emission spectrum 
at 10 7 or 10 6 - 3 K(see text). Incident and H/3 flux in erg cm s . Equivalent widths relative to nebular continuum only. 
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Fig. 1. — Initial configuration (t = 0) of simulated accretion disk before recoil of the central black hole (see text). Inner and outer radii 
are at 10 and 50 in the dimensionless units of the model. Disk rotates counterclockwise seen from above. 
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Fig. 2. — Test particle simulation of accretion disk at time t = 1 after recoil of the central black hole (see text). In the dimensionlcss 
units of the model, the recoil has = 65 and is inclined 45 degrees to the initial disk axis. This corresponds to 65% of the innermost 

orbital velocity. The pre-kick Keplerian orbital period is P b = 2.28 at the radius R b = 23.6 where i>K e p = ''kick- The orbital period at the 
inner and outer boundaries of the initial disk are P = 0.63 and 7.0, respectively. 
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Fig. 3. — Simulated accretion disk at time t = 2.5 = l.lPb after recoil of central black hole (see text). 
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Fig. 4. — Predicted light curves for recoil flare (see text). Black curve with sustained luminosity is power per unit time for particle 
collisions in numerical simulation; power at later times is overestimated because of repeat collisions. Green curve with declining power is 
based on post-kick orbital period P of each bound particle, assuming that an energy (l/2)r?ipD^ ick is dissipated at t = P. Smooth red curve 
is analytic solution for axial kick. Each curve is normalized to Pb for the time axis and Lo = (l/2)MbiijJ ick Pb for the luminosity axis. For 
reference, P b = 10 3 4 yr and Lo = 10 45 7 erg s" 1 for M BH = 10 8 M Q , M = 10~ o s Mq yr -1 and ii kick = 1000 km s -1 (Equations g] and 
0. 
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Fig. 5. — Spectrum of optically thin gas in coronal equilibrium at T = 10 7 K in wavelength range 5 to 20 A (see text). Note strong 
emission lines of Fe XX and neighboring ions. Abscissa is wavelength in A in bins 1500 km s _1 wide; ordinate is \F\ in arbitrary units. 
The intrinsic line width for the plot is set at 1000 km s _1 . 
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Fig. 6. — Spectrum of photoionized precursor for Tg^ock = 10 7 K in wavelength range 800 to 2000 A (see text). Note strong emission 
lines of O VI, Lya, and C IV. Flux is \F\ in arbitrary units. 
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Fig. 7. — Spectrum of photoionized shock precursor for Tg^ock = 10 7 K in wavelength range 2600 to 5000 A (see text). Flux is \F\ in 
arbitrary units. Note lines of H/3, He II A4686, and Mg II, and Balmer continuum in emission. 



